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LOW-AMPLITUDE DAMPING-IN-PITCH CHARACTERISTICS OF FOUR 
TATLLESS SWEPT WING-BODY COMBINATIONS AT MACH 
NUMBERS FROM 0.85 TO 1.30 AS OBTAINED 
WITH ROCKETT-POWERED MODELS 


By Charles T. D'Aiutolo 
SUMMARY 


A free-flight investigation was conducted to determine the damping- 
in-pitch characteristics of four rocket-powered wing-body combinations 
between the Mach numbers of 0.85 and 1.50. А11 models were tested at the 
same center-of-gravity location (17 percent mean aerodynamic chord) and 
had identical bodies thereby allowing a direct comparison of the aerody- 
namic characteristics of the various wings. The wings tested were an 
aspect-ratio-h, 45° swept composite-plan-form wing (derived from an 
aspect-ratio-6 wing by increasing the area behind the trailing edge of 
the inboard 40 percent of the wing), an aspect-ratio-h, 45° swept highly 
tapered wing, an aspect-ratio-3, 52.59 swept highly tapered wing, and an 
aspect-ratio-3, 609 swept highly tapered wing. Тһе composite wing had 
NACA 65A006 airfoil sections parallel to the free stream, whereas the 


swept highly tapered wings had 65A004 airfoil sections parallel to the 
free stream. 


The results indicated that, although all models were statically stable 
longitudinally throughout the Mach number range investigated, the dynamic 
stability of the models was low and differed appreciably for each model. 
The total damping factor and the slope of the lift curve were about the 
same order of magnitude resulting in values of the rotational damping-in- 
pitch derivatives that were either negative (stable) or positive 
(unstable). 


A comparison of the damping was made between a delta wing in combina- 
tion with a body and a swept wing in combination with a body and the 
results indicated that at M = 1.16 the delta-wing configuration had better 
damping characteristics. ы 


ы йрек у m S À 
L ta c 23 


— a -— n ee 


2 | у. 2 NACA RM 15%110 


INTRODUCTION 


Recently & series of systematic experimental investigations has been 
conducted on tailless aircraft configurations to verify the theoretical 
predictions of low dynamic longitudinal stability and regions of dynamic 
instability of these designs (see, for example, refs. 1 to 5). Inasmuch 
as most of this work was concerned with the delta-wing configurations and 
very little information is available concerning swept-wing configurations, 
this investigation was undertaken to determine experimentally the regions 
of dynamic instability at transonic speeds and low supersonic speeds of 
some tailless swept-wing airplane configurations. 


This report is a continuation of the investigation reported in ref- 
erence 5 and contains the results from the flight tests, conducted by the 
Pilotless Aircraft Research Division, of four rocket-powered wing-body 
combinations. Two of the wings used in these configurations had 45° of 
sweep and aspéct ratios of 4, the only difference being that one was of 
composite plan form (resulting from & search for an improved design suit- 
able for high-speed flight, see ref. 6) and the other was highly tapered. 
The third and fourth wings were both of aspect ratio 5 and highly tapered, 
with one having 52.59 of sweep and the.other 60°. ‘The composite wing had 
NACA 65A006 airfoil sections parallel to the free stream, whereas the 
swept highly tapered wing had NACA 65A004 airfoil sections parallel to 
the free stream. The data are presented over a Mach number range of about 
0.85 to 1.30 corresponding to a Reynolds number range (based on respective 
mean aerodynamic chords) of about 7 x 106 to 14 x 106, respectively. 


The models were flown at the Langley Pilotless Aircraft Research 
Station at Wallops Island, Va. The static and dynamic longitudinal 
stability characteristics of the models were determined by analyzing 
the oscillations produced in pitch by firing small pulse rocket motors 
which were mounted to provide thrust normal to the longitudinal axis of 
the models. The drag characteristics of the models and the stability 
characteristics were determined from the deceleration portion of the 
flights. 


SYMBOLS 


V velocity of flight, ft/sec 


V 


M Mach number, | ————————— 
Speed of sound 


S total wing area, sq ft 
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Cl 
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mean aerodynamic chord, ft 
aspect ratio 


cross-sectional area of any longitudinal station of complete 
model, sq ft 


body length, ft 
distance along body measured from nose, ft 
&ir density, slugs/cu ft 


Reynolds number (based on mean aerodynamic chord of respective 
models) 


free-stream dynamic pressure, 1b/sq ft 

total damping factor (logarithmic decrement of pitch oscillations) 

period of short-period oscillations, sec 

reduced-frequency parameter (based on respective mean aerodynamic 
chords of models), © 


2V 


lift coefficient, Cy cos a - Cc sin a 


À 
normal-force coefficient, Ax W/S 
Б а 


acceleration normal to reference axis as obtained from 
accelerometer 


W/S 


chord-force coefficient, - r: E 


acceleration along reference axis as obtained from accelerometer 
positive forward 
model weight, lb 
Ò 
slope of lift curve per degree, zn 


pitching-moment coefficient, — 


ү е 
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Ст, static-stability derivative per degree, Lu 
ӘС 
Сша = qE) per radian 
2V 
Сва = Sa per radian 
oí С. 
2V 
Cp drag coefficient, Cc cos a + Cy sina 
a angle of attack, measured from flight path to fuselage reference 
line, deg 
а = —— 202 radians/sec 
97.35 dt 
Ө angle of pitch, measured from horizontal to fuselage reference 
line, radians 
6 = I, raaians/sec 
dt 
t time, sec 
w frequency of short-period oscillation, radians/sec 
= —®— relative-density factor 
рве 
Iy moment of inertia in pitch, slug-ft? 
m mass, slugs 


MODELS AND APPARATUS 


Model Description 


The general arrangements of the models are shown in figure 1, and 
a photograph of one of the models is presented in figure 2. Geometric 
characteristics of the models are presented in table I. Model 1 had an 
aspect-ratio-h, 45° swept composite plan-form wing (derived from an aspect- 
ratio-6 wing by. increasing the area behind the trailing edge of the inboard 
ҺО percent of the wing) with a taper ratio which varied from about 0.7 for 
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the outboard panel to about 0.5 for the inboard panel and incorporating 
NACA 65A006 airfoil sections parallel to the free stream. Model 2 had 
an aspect-ratio-h, 45° swept highly tapered wing, model 3 had an aspect- 
ratio-3, 52.59 swept highly tapered wing, and model ! had an aspect- 
r&tio-5, 60° swept highly tapered wing. The wings of models 2, 3, and 4 
incorporated NACA 65A004 airfoil sections parallel to the free stream. 
Note that the wings of models 1 and 2 were identical in aspect ratio and 
sweep, but differed in plan-form shape and section thickness, whereas the 
wings of models 3 and 4 were identical except for wing sweep. 


Each model consisted of a basic fuselage to which was attached the 
wing to be tested. The fuselage was essentially a body of revolution 
with rather large wing-root fairings for structural purposes and consisted 
of an ogival nose section which contained the telemeter and a cylindrical 
body section which contained the wing mount, necessary fairings, the 
vertical tail, and the sustainer rocket motor. Construction of the fuse- 
lage was of aluminum alloy with magnesium skin. 


The wings of model 1 were constructed of wood with sheet-steel inlays 
(for structural purposes), whereas the wings of models 2, 3, and 4 were 
constructed of solid aluminum. All wings were mounted on the fuselage 
(as shown in fig. 1) with the resultant center of gravity located at 
17 percent of the mean aerodynamic chord during the deceleration portions 
of the flights. 


Each model contained a Cordite sustainer rocket motor and was boosted 
by a lightweight 5-inch HVAR rocket’ motor. А1] models with their 
boosters were launched from a rail-type launcher (as shown in fig. 3) at 
an angle of approximately 5°. 


The cross-sectional-area distribution of each of the models is 
presented in figure 4, for possible correlation of drag results. 


Instrumentation 


Each model contained a standard 4-channel NACA telemeter which trans- 
mitted continuous flight measurements of angle of attack (measured by a 
vane-type instrument located on a sting forward of the nose of the 
models), normal acceleration of the center of gravity, longitudinal accel- 
eration, and pitot stagnation pressure (measured by a tube located on a 
strut below the fuselage of the models). (See fig. 1.) 


The position of the models in space was determined from an SCR 58h 
tracking radar and the velocity of the models was obtained by use of the 
Doppler velocimeter radar. Atmospheric data were obtained from a radio- 
sonde released just prior to each of the model flights. Fixed and 
tracking motion-picture cameras were used to observe the conditions of 
the models during the major portion of the flights. 


amt 
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TEST AND ANALYSIS 


Test 


The data for each model were obtained during the decelerating portion 
of the flights. The damping-in-pitch data were determined by disturbing 
each of the models in pitch by & series of four small rocket motors pro- 
viding thrust normal to the longitudinal axis of the model and located 
near the rear of the model. These rocket motors were timed to fire in 
sequence during the decelerating portion of ihe flights; however, the 
first small rocket motor to be fired in each of the four models fired 
prematurely during the latter part of sustainer motor burning. Also, 
one of the remaining three small rocket motors in model ! failed to fire, 
so that damping data from the flight test of model 4 were obtained from 
the firing of only two of these small rocket motors. Thus, in most 
instances, reliable damping data could be only obtajned from three of 
the disturbances. 


Time histories of angle of attack, lift coefficient, and Mach number 
covering the decelerating portions of the flight and the times at which 
the small pulse rockets were fired for each of the four models are shown 
in figure 5. Also shown are the envelopes drawn for each of the oscil- 
lations in pitch that were caused by the firing of the small rocket motors. 
The static and dynamic longitudinal stability derivatives were obtained 
from these oscillations. 


Data, obtained from spinsonde records, indicated that for each of 
the four models the rate of roll was approximately zero. 


Figures 6, 7, and 8 present, respectively, the variation of air 
density, velocity, amd dynamic pressure with Mach number for each of the 
tests. These quantities are presented so that a possible correlation of 
the data obtained from these tests with data obtained from other tests 
may be made. The range of Reynolds numbers of the tests is plotted against 
Mach number in figure 9 where the Reynolds number is based on the respec- 
tive mean &erodynamic chords of the models. 


Accuracy 
The maximum probable errors in the basic coefficients and angle of 
attack are as follows: 
M = 0.85 M = 1.25 
Dn ә oer LUE ce OR 0н ecc ж Q ¿g e WZ жок жок 0.007 0.005 


ACO: won GR ж-ш эй уй. EDO: жошоо Ga. CE. SEL dE c) od dee ж 0.007T 0.005 
Ad, deg. cxx ж $09 Wo; ж 990€. ie o e cw ДА M des 16 0.20 о .20 
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These errors are largely due to errors in instrumentation which is a 
function of the full-scale range of the instruments and errors due to 
flight conditions. The Mach number as determined from Doppler radar 
velocity amd radiosonde measurements is accurate to better than 2 percent. 


Analysis 


The analysis of the motions of the models presented in this report 
are based on the assumption of constant coefficients in the differential 
equations of motion. (See ref. T.) 


Ihe slope of the lift curve Cry, was found by plotting CI, against 


a as obtained from the flight records during an oscillation and graphi- 
cally measuring the slope. The Cm, data that are presented in this 
paper were calculated by use of the expression 


we Salt : 


which is based on a single-degree-of-freedom analysis instead of the more 
usual expression 


-— 1 (85) + >| (2) 
27.5998 V Р 
which is based оп а two-degree-of-freedom analysis. For tailless соп- 
figurations, the contribution of the total damping factor b to Cng, 
is negligible when compared to the frequency contribution, so that the 
single-degree-of-freedom expression given (eq. (1)) allows an accurate 
determination of Cm,. Values of the total damping factor b меге 
determined by the method presented in the "Analysis" section of refer- 
ence 5, while values of Cmq + Cm& were determined by the method pre- 
sented in appendix A of reference T. 


The maximum value of the reduced-frequency parameter k = == (based 


on respective mean aerodynamic chords) that was determined from these 
tests was O.025. Since this value is small, it is believed that the 
effects of the frequency of the oscillations of the models on the damping 
in pitch are not important im the determination of the damping-im-pitch 
derivatives amd that the method of reducing the data used in this paper 
gives good results for these derivatives.  Unpublished calculations, 
comparing the damping in pitch computed for terms to the order of К 
with the damping in pitch computed for terms to the order of k, indicate 
that for k < 0.025 the difference is about 1 percent. These results 
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indicated that terms of higher-order frequency are not required in the 
estimation of the damping-stability derivatives. 


RESULTS AND DISCUSSION 


The stability parameters of the models presented in this paper were 
determined from the coasting phase of the flights. Each of the models 
was tested with the center of gravity located at 17 percent of the mean 
aerodynamic chord behind the leading edge of the mean aerodynamic chord. 


Trin 


The trim characteristics of each model are shown in figure 5. 
Model 3 (A = 3, 52.5° sweep) experienced a small trim change in the tran- 
sonic speed range; however, it was the largest trim experienced by any of 
the models tested. At high subsonic speeds, model 3 flew at a trim angle 
of attack of approximately 1.0° corresponding to a trim lift coefficient 
of approximately 0.02 and experienced a nose-down trim change at transonic 
speeds. At supersonic speeds, model 5 flew at a trim angle of attack 
of approximately 0.5° corresponding to a zero-lift condition. 


Lift 


The variation of the slope of the lift curve for the four models is 
presented in figure 10. Тһе test points represent data that were deter- 
mined from the oscillations produced by the firing of the small pulse 
rockets and, inasmuch as limited data were obtained, no attempt was made 
to fair curves through the test data. 


Models 1 and 2 (identical aspect. ratio and sweep, but different in 
plan-form shape and section thickness) show about the same variation with 
Mach number; although, at transonic speeds and supersonic speeds, Clay 
of model 2 is somewhat lower than the CI, of model l. These differences 
are not due to wing flexibility since preflight static wing twist indi- 
cated that models 1 and 2 had about the same degree of flexibility at 
transonic and supersonic speeds. 


The effect of sweep on Cla may be seen by comparing model 5 
(52.59 sweep) with model 4 (60° sweep) since, except for wing sweepback, 
the models were the same. Also, both models had about the same degree 
of wing flexibility. The Сї data of model h are somewhat lower than 
the CI, data of model 5 as would be expected since the wing sweep of 
model 4 was greater than the wing sweep of model 3. 
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Static Longitudinal Stability 


The static longitudinal stability for each of the four models is 
presented in figure ll where it is seen that Cm, increases with 
increasing Mach number through the transonic speed range, then decreases 
somewhat as the Mach number becomes supersonic. From a comparison of 
models 1 and 2 (same aspect ratio and sweep, but different in plan-form 
shape and section thickness), it is seen that at high subsonic and tran- 
sonic speeds model 1 has greater static longitudinal stability than 
model 2, while at supersonic speeds the reverse is true. A comparison of 
model 3 with model 4 (same aspect ratio, but different sweepback) shows 
that model 4 has slightly greater static longitudinal stability at high 
subsonic and transonic speeds than model 4, but at supersonic speeds the 
reverse is true. Note that the static stability of models 2 and lh 
increase to a higher Mach number before decreasing than models l and 3. 


Since all wings were mounted on identical bodies and were tested at 
the same center-of-gravity location with respect to the mean aerodynamic 
chord, a comparison of the aerodynamic-center location of the various 
wings can be made. These values of the aerodynamic-center location appear 
in the following таре. 


Aerodynamic-center location in percent 
mean aerodynamic chord for - 


Subsonic . . 
Transonic . 





Note that generally for all models as the Mach number increases the 
&erodynamic-center location moves rearward; however, model 3 indicates 

a slight forward movement as the Mach number increases from high subsonic 
speeds to transonic speeds. It is noted that this forward movement of 
the aerodynamic-center location occurred while model 3 experienced a 
nose-down trim change (see fig. 5(c)). 


Dynamic Longitudinal Stability 


Total damping factor.- The total damping factor b was obtained 
from the time histories of the pitch oscillations (see fig. 5). This 


total damping factor includes the contribution of moment due to motion 
along a curved path at constant angle of attack Cg, the moment due to 


plunging motion (vertical acceleration) Сых, and the translational 
GONE THENT IAS, 
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effect of Ст. The variation of b with Mach number for each of the 


four models is presented im figure 12. By use of the method presented 
in the "Analysis" section of reference 5 it was possible to determine 
numerous instantaneous values of b at high subsonic speeds for model 1 
and at low supersonic speeds for models 2, 5, and h such that the solid 
curves presented in figure 12 represent values that were determined from 
experimental test data points only. The dashed curves represent fairings 
between other data points that were obtained for models 1, 2, and 3. 


Figure 12 shows that plan-form shape has a pronounced effect on the 
total damping factor for the low amplitudes reported in this paper. A 
comparison of models 1 and 2 (identical aspect ratio and sweep, but dif- 
ferent in plan-form shape and section thickness) shows that for the 
composite-plan-form configuration (model 1) the total damping factor Ъ 
decreases abruptly then increases abruptly for small increments of Mach 
number at high subsonic speeds, while for model 2 this abrupt decrease 
followed by an abrupt increase in b over a small range of Mach numbers 
occurs at low supersonic speeds. These abrupt changes were determined 
from amplitudes that were less than 0. 5 1 in angle of attack and therefore 
may not be significant. 


The total-demping-factor data for model 5 (52.59 sweep) and model hà 
(609 sweep) are also shown in figure 12. Inasmuch as limited data were 
obtained, it was not possible to determine the effect of sweep on ihe 
total damping factor. 


Rotational damping-in-pitch derivatives.- The rotational damping-in- 
pitch derivative for each of the four models was determined by the fol- 


‘lowing expression: 


_ EVIY(/ gb , Clg 
Сша + Cm& = (57.598 + E) (5) 


where m is the mass of the model in slugs. It is seen from the fore- 
going equation (5) that the determination of the value of the rotational 
damping-in-pitch derivatives is dependent upon differences in the total 
demping factor b and the translational effect due to Ст. Im the model 


tests reported in this paper (tailless configurations) the total damping 
factor and the translational effect of CI, are about the same order of 
magnitude resulting in small values of Cmq + Ст that are either nega- 
tive (stable) or positive (unstable). (See figure 13.) 


Effect of plan-form shape on the rotational damping-in-pitch 
derivatives.- A comparison is made in figure 14 between an aspect-ratio-3, 


tailless delta wing in combination with a body and an aspect-ratio-5, 
52.5° swept highly tapered wing in combination with a body (model 3) to 


quum» 
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determine the effect of plan-form shape on the rotational damping-in-pitch 
derivatives. The data for the delta-wing configuration were obtained from 
reference 5. 


From figure 14, it appears that there is little effect of plan-form 
shape on the rotational damping-in-pitch derivatives at M = 0.915 and 
M = 1.05; however, at М = 1.16 the delta-wing configuration has better 
damping characteristics than the swept-wing configuration. Reference 1 
also shows that the damping in pitch of the delta-wing plan form is 
superior to the swept-wing plan form at supersonic speeds. 


Drag 


The variation of the drag coefficient with Mach number for each of 
the four models is presented in figure 15. The effect of wing sweep may 
be seen by comparing models 3 and 4+. Model 4 has lower drag throughout 
the Mach number range than model 3 as would be expected since model 4 has 
a greater wing sweep than model 5. An indication of the effect of plan 
form may be made by comparing models 1 and 2 where it is seen that the 
composite plan form has considerably better drag characteristics than the 
conventional swept plan form even though the wing of the composite plan 
form was somewhat thicker than that of the conventional swept plan form. 


CONCLUSIONS 


From the results of the free-flight tests of four wing-body combina- 
tions incorporating identical bodies and tested at the same center-of- 
gravity location of 17 percent mean aerodynamic chord and consisting of 
an aspect-ratio-4, 45° swept composite wing, an aspect-ra&atio-h, 45° swept 
highly tapered wing, an aspect-ratio-3, 52.5° swept highly tapered wing, 
and an aspect-ratio-3, 60° swept highly tapered wing, the following 
conclusions may be stated: 


1. All models were statically stable throughout the Mach number range 
investigated (M = 0.85 to 1.30). 


2. For the low amplitudes reported in this paper, the total damping 
factor was affected by plan-form shape, whereas the rotational damping- 
in-pitch derivatives of the models were low throughout the Mach number 
range investigated and differed appreciably between models. Inasmuch as 
the damping-in-pitch derivatives were low, the total damping factor con- 
sisted mainly of the contribution of the slope of the lift curve. 


| 
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3. At М = 1.16 and for the wings with aspect ratio 5, the delta- 


wing configuration had better rotational damping-in-pitch characteristics 
than the swept-wing configuration. 


Langley Aeronautical Laboratory, 


National Advisory Committee for Aeronautics, 
Langley Field, Va., August 27, 1954. 
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TABLE I.- GEOMETRIC CHARACTERISTICS OF THE MODELS 


Wing: 
Total area, sq ft. 
Span, ft š : 
Aspect ratio: 


Mean aerodynamic chora, Tb dex» ux u 
Sweepback of querter-chord, deg. . 


Dihedral, deg : 

Taper ratio - . 

NACA airfoil sections parallel 
to free stream 


Fuselage: 


Miscellaneous: 
Model weight, lb 
Moment of inertia in pitch, 
Iy, slug-ft2 ... 
Center-of-gravity position, 
percent M.A.C. . 


Wing loading, lb/sq ft ...... 


Relative density factor, u - 


5.76 
l. 82 
М 
1.59 
L5 


О 
80.715, bo.288 


65A006 


12.00 
10.70 


128.25 


Model 3 


Model 4 
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Figure l.- 


Pitot stagnation tube 


15.00 


(а) Model 1. А = 43 airfoil sections, NACA 654006. 


General arrangement of models. All dimensions in inches unless 
ctherwise noted. 
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(b) Model 2, A= N; airfoil sections, NACA 65AO0h. 


Figure 1.- Continued. 
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А = 3; airfoil sections, NACA 65A004. 


Figure 1.- Continued. 
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Figure 1.- Concluded. 


(a) Model h, А = 5; Airfoil sections, NACA 65A00h. 
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Figure 2.- Photograph of one of the models. 
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Figure 5.- Photograph of one of the models on the launcher. 
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Figure 4.- Longitudinal distribution of cross-sectional area. 
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(a) Model 1. 


Figure 5.- Time history of flights showing envelopes of the oscillations, 
trim values, and the times where the pulse rockets were fired. 
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(b) Model 2. 


Figure 5.- Continued. 
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(c) Model 3. 


Figure 5.- Continued. 
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Figure 5.- Concluded. 
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Figure 6.- Variation of air density with Mach number. 
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Figure 7.- Variation of velocity with Mach number. 


NACA RM 154110 СОЛЛА, ` 27 


Model 2 
1600 A=), Z 
NACA 65A00l 


T Ail | | 
1b/sq a TT TT TAR з 
HHHH 945 LE 
1200 HH PH 
ti | A a 
1000 КЕГЕ EAE чини 





Figure 8.- Variation of dynamic pressure with Mach number. 
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with Mach number. 


dynamic chord 


Figure 9.- Variation of Reynolds number, based on respective mean aero- 
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(b) Models 3 and }. 


Figure 10.- Variation of the slope of the lift curve with Mach number. 
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(b) Models 3 and L. 


Figure ll.- Variation of the slope of the pitching-moment curve with Mach 
number. 
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(b) Models 3 and h. 


Figure 12.- Variation of the total damping factor with Mach number. 
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(b) Models 5 and h. 


Figure 15.- Variation of the rotational damping-in-pitch derivatives with 
Mach number. 
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Figure 1}.~ Effect of plan-form shape on the rotational damping-in-pitch 
derivatives. 
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Figure 15.- Variation of the drag coefficient with Mach number. 
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